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SUMMARY 


Goals and objectives 


The long range goal of this program is to develop an improved understanding of 
phenomena of importance to directional solidification, to enable explanation and 
prediction of differences in behavior between solidification on Earth and in space. 
Current emphasis is on determining the influence of perturbations on directional 
solidification. 


Progress 

Mohsen Banan completed his Ph.D. thesis and took a job with MEMC in St. Peters, 
Missouri, to work on Czochralski growth of silicon. Excerpts from his thesis are 
included here as Appendix A. A new graduate student, Mr. Weijun Yuan, began work 
in January. Until August he is being supported by the Department of Chemical 
Engineering as a teaching assistant. He has been learning techniques for preparing 
ampoules, directional solidification, and metallography. 

The paper "Experimental observation of the influence of furnace temperature profile on 
convection and segregation in the vertical Bridgman crystal growth technique," 
presented at the Dresden LAF meeting last October, was accepted for publication in Acta 
Astronautica. A paper on the influence of vibration on the microstructure of lead-tin 
eutectic was accepted for publication by the Journal of Crystal Growth. It is reproduced 
here as Appendix B. 

Mohsen Banan directionally solidified InSb-GaSb alloys while the ampoule was vibrated 
or electric current pulses were applied through the ingot. In all cases the axial 
composition profile corresponded closely to that expected with complete mixing in the 
melt and equilibrium at the freezing interface. Interface breakdown occurred near the 
end of the ingot, as indicated by a finer microstructure, large fluctuations in 
composition, and microcracks. Application of current pulses caused the onset of these 
conditions much sooner, indicating that these pulses increase constitutional 
supercooling. Vibration appeared to decrease significantly the number of grains and 
increase slightly the number of twins, which nucleated at the ampoule walls. 

James Zhou directionally solidified Te-doped InSb with periodic rotation of the ampoule 
about its axis (ACRT). Striations were produced showing that the growth rate varied 
significantly down the ingot, the interface became facetted near the end, the freezing rate 
was different when rotation was on than when it was off, and the freezing rate varied as 
the crucible rotated. 


i 


Plans 


A new graduate student, Weijun Yuan, will follow up on Banan's observations. He will 
determine the reproducibility of the influence of vibration on microstructure and, if the 
effect is real, try to determine the mechanism causing it. James Zhou will perform 
solidification experiments with both ACRT and current interface demarcation in order to 
determine the variation of growth rate and interface shaped during each ACRT cycle. 


Publication 


R. Caram, M. Banan and W.R. Wilcox, "Directional solidification of Pb-Sn eutectic with 
vibration," J. Crystal Growth (in press). (See Appendix B.) 
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I. INFLUENCE OF IMPOSED PERTURBATIONS ON DIRECTIONAL 
SOLIDIFICATION OF InSb-GaSb ALLOY SEMICONDUCTOR 

Mohsen Banan and Weijun Yuan 


Summary 


The objective of this project is to determine the influence of imposed perturbations on 
compositional homogeneity and microstructure of directionally solidified alloy 
semiconductors. Solid solutions of InSb-GaSb were solidified in a vertical Bridgman - 
Stockbarger apparatus. Two different approaches were undertaken to impose 
perturbations in the growth system: 

1. Periodic reversed of the direction of electric current passed axially through the 
melt-solid configuration. 

2. Axial vibration of the growth ampoule at 10-100 FIz frequency and 0.05-1.0 mm 
amplitude. 

In x Gai-xSb ingots with x = 0.2 initial feed composition were directionally solidified 
using a vertical Bridgman-Stockbarger technique at 8 mm/day and 25-30 C/cm axial 
temperature gradient in the empty furnace. Energy dispersive x-ray spectrometry was 
used for compositional analysis of the ingots. The microstructure was examined using 
optical and scanning electron microscopy. 

Severed ingots were grown with axial vibration of the ampoule. Vibration at 10-100 Hz 
and 0.05 to 1.0 mm amplitude was applied to the growth ampoule during growth. 
Several ingots grown with vibration had fewer grains and more twins than one solidified 
without vibration. For example, ingots grown under vibration at 10 Hz of 1.0 mm 
amplitude and 20 Hz of 0.5 mm amplitude consisted of 2 grains compared to 7 to 9 
grains in the ingot solidified without vibration. 

Both with and without vibration the axial composition profiles corresponded to those 
expected for complete-mixing in the melt. The composition across the sample (radially) 
was uniform. 

An ingot was grown with 15 amp/ cm ^ of electric current passed axially through the 
solid and the melt. The current was alternatively changed in direction every 25 sec 
during growth. The current-induced perturbations caused an earlier transition to a fine 
grain structure with large fluctuations in composition and micro-cracking, all indicative 
of interface breakdown due to constitutional supercooling. 

Mr. Banan completed his Ph.D. thesis, portions of which are in Appendix A. A new 
student, Mr. Yuan, will determine the reproducibility of the influence of vibration on the 
microstructure. 
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II. USE OF CURRENT INTERFACE DEMARCATION DURING 
DIRECTIONAL SOLIDIFICATION OF Te-DOPED InSb WITH ACRT 

James Zhou 


Summary 


The objective of this project is to use the Current Interface Demarcation (CID) technique 
during the solidification of Te-doped InSb to determine the influence of Accelerated 
Crucible Rotation Technique (ACRT) on the growth. We expect to create several CID 
striations during each ACRT cycle by passing current pulses. Observing these striations 
will allow us to determine the change of growth rate and interface shape during the 
ACRT cycle. 

Tellurium-doped InSb was solidified with ACRT at different rotation rates and with 
different dopant concentrations, but without CID. Thick striations were induced by 
turning on and off the rotation. In addition, there was a thin rotational striation for each 
revolution of the ampoule within the sections solidified with the rotation on. The 
rotation rate and dopant concentration influenced the appearance of the striations. 

The growth rate was calculated by measuring the distance between adjacent striations. 
It decreased rapidly at the beginning of the ingots, became nearly constant, increased 
slightly and then diminished to near zero at the end. The striations changed direction 
when they crossed grain and twin boundaries. The interface was totally faceted at the 
tail end of the ingots. Under the microscope, the striations caused by spin-up looked 
different from those caused by spin-down. The length solidified during spin-up was 
also different from that during spin-down, suggesting that the growth rate was different 
during spin-up and spin-down. 

A novel ampoule was made to allow passage of current pulses during solidification with 
ACRT. Solidification will be performed with both ACRT and CID. 

An MS thesis will be finished during the next reporting period. 


II.A. Introduction 


Developed by Scheel et al. [1] in the early 1970s, the Accelerated Crucible Rotation 
Technique (ACRT) has been employed in various crystal growth processes. For 
example, with ACRT crystals grew 1000 times larger than without ACRT in the 
high-temperature flux growth of GdAlOg [2]. Vertical Bridgman-Stockbarger (VBS) 
growth of Cd x Hg|. x Te with ACRT was studied extensively by Capper et al. [3-6]. 
Improved compositional homogeneity and larger grains were obtained. Demarcation of 
the interface by quenching showed that ACRT caused the growth interface to become 
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less concave [6]. 


Gray [7] at Clarkson has been studying the influence of ACRT on InSb-GaSb 
solidification. Both with and without ACRT, his axial composition profiles always 
corresponded to complete mixing in the melt and equilibrium at the interface. ACRT 
appears to have had only a slight influence on the microstructure. 

The mechanism for the above behavior during solidification with ACRT is not clear. 
Periodic acceleration and deceleration of crucible rotation leads to effective stirring of the 
melt. Thus improved cross-sectional compositional homogeneity is expected in the 
solidified ingots. Several mechanisms have been proposed to explain grain 
enhancement by ACRT, including a more favorable interface shape, backmelting and 
regrowth, and "scrubbing" of nuclei. 

Many crystal growers have suggested that the periodic nature of ACRT may produce a 
periodic variation in the mixing of the melt and a periodic variation of the temperature 
adjacent to the interface. Both of these would be expected to cause compositional 
striations to be grown into the ingots. Larrousse [8] at Clarkson has confirmed the 
periodic nature of mixing resulting from spin-up and spin-down by limiting current 
density electrochemical experiments. From measurements of temperature distributions 
during high temperature solution growth, Tolksdorf and Welz [9] concluded that there 
were temperature oscillations in synchronism with the program of variable rotation. 

These observations imply that the composition of the grown solid would also fluctuate, 
producing striations. These striations might be useful as time markers and to reveal the 
interface shape. Striations were found when Gornert et al. [10] investigated the flux 
growth of YIG with ACRT. In other systems, however, the high temperature present in 
the solid near the interface could cause these striations to diffuse away before the solid 
reaches room temperature, and so never be seen. Gray et al. [11] developed a one 
dimensional, time-dependent model for the diffusional decay of striations. Neither 
Capper et al. nor Gray found striations in their systems. 

For our Te-doped InSb system with ACRT of 10 s on and 10 s, we estimate 
D/V^t c = 2.75 x 10-5, w here D is the diffusivity of Te in solid InSb at its melting point, V 
is the solidification rate, and t f is the ACRT cycle time. Using the theory of Gray et al 
[11] we estimate the distance for 99.9% decay to be 35 cm, which means that striations 
caused by ACRT should not diffuse away in the present experiments. 

Passage of current through the melt-solid interface causes Peltier cooling or heating, 
depending on the direction of the applied current [12,13]. Since tellurium is rejected by 
the growing solid, an increase in growth rate should lead to a Te-rich striation in the 
solidified material. Generation of several striations within the material grown in one 
ACRT cycle would enable us to determine the variation of growth rate and interface 
shape during that cycle. 


3 



II.B. Experimental Methods and Results 


Figures 1 and 2 are schematic diagrams of the ampoule and the furnace system that 
allow passage of current pulses during ACRT. The ampoule is about 85 cm long, 9mm 
ID and 11mm OD. The current pulses pass from a graphite electrode at the bottom, a 
molybdenum lead, another graphite electrode, through the ingot, another molybdenum 
lead, and finally to the graphite electrode at the top of the ampoule. Both the top and 
bottom graphite electrodes are contacted with copper brush electrodes which attach to 
the current pulse generator. A bearing is mounted on each end of the furnace to align 
the ampoule. 

Several growth runs with ACRT only and one with both ACRT and CID were 
performed. The growth conditions are listed in Table 1. 

For the first run Rl, an ACRT rotation rate of 20 rpm was applied at the beginning. After 
about 10 hours, the ampoule was lowered down 1 mm and the rotation rate was 
changed to 80 rpm. The cycle time was 10 seconds with rotation on and 3 minutes with 
it off. The resulting ingot was sliced longitudinally, mechanically polished, and 
chemically etched. This section was observed under the microscope at magnifications 
from about 40x to 400x. There were no striations at the beginning portion. The mark 
caused by the fast lowering was about 3.1 mm from the front of the ingot. After the 
mark, striations were visible to the naked eye when the sample was held at a certain 
angle, which depended on the orientation of the grains and the light source. These same 
striations were not observable under brightfield, darkfield, and even Nomarski 
microscopy. Striations became more pronounced at the tail end of the ingot. Figures 3 
and 4 show the striations in the tail end of ingot Rl. These striations showed that the 
melt-solid interface was not flat at the end of the growth. It changed its direction 
sharply across boundaries. The dark region in Figure 5 might be oxidized material. 
Figure 6 is a scanning electron micrograph showing striations in grains about 1.3 mm 
from the end. The distance between the striations showed the growth rate decreased as 
the end was approached. 

In order to learn the reason for the absence of striations at the beginning of Rl, ingots R2 
and R6 were grown under identical conditions, except an 80 rpm rotation rate was 
applied for the entire time of these runs. The Te concentration in R6 (5 x 10 iy / ott) was 
10 times that in Rl and R2. There were still no striations at the beginning of ingot R2. 
However, the first striation appeared much earlier, 9 mm from the beginning, as shown 
in Figure 7. In ingot R6 striations were observed from the very beginning of the ingot, 
indicating that a larger tellurium concentration is required to reveal striations under the 
conditions employed. 

There were two dominant striations for each ACRT on-off cycle, corresponding to the 
start and stop of rotation. One interesting observation, as shown in Figure 8, was that 
there were fine striations produced during the 10 seconds the rotation was on. One such 
stration was produced for each revolution of the ampoule, indicating that the thermal 
environment was sufficiently non-uniform to cause a periodic variation in freezing rate 
at each point. 

In runs R3 and R4 ACRT was used with the time of roation equal to the time without 
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rotation. The striations caused by spin-up looked different from those caused by 
spin-down, as shown in Figure 9. The length solidified during spin-up was also 
different from that during spin-down, indicating that the growth rate was different 
during spin-up and spin-down. 

Figure 10 shows coupled striations in ingot R1 that indicate that the average growth rate 
during rotation was twice that without rotation. Further investigation will be 
undertaken. 

One growth run was performed with both ACRT and CID. The resulting ingot is being 
studied. 


II.C. Future Work 

• More growth runs with both ACRT and CID. 

• Several growth runs with ACRT only to ascertain whether turning on the rotation 
causes a larger or smaller freezing rate than turning off the rotation. 

• Some experiments with CID only to investigate the influence of the length of the 
insulation zone and the length of ingot on the growth. 
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Table 1: Growth Conditions 
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Figure 2 EXPERIMENTAL SETUP 
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Figure 3: Ingot R1 with high magnification shows Stria- 

tions induced by ACRT. (250x) 
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Figure 4: Striations changed direction when they crossed 

the boundaries. 1.3 mm from tail end of ingot R1 (225x) 
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Figure 5: Faceted tail end of ingot Rl. Dark region might 

be oxidized material. Growth direction was shown by arrow. 
(40x) 
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Figure 6: Striations at the end of ingot R1 showed the 

growth rate diminish to near zero when the growth finished 
(lOOx) 
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Figure 8: Rotational striations about 10 mm from begin 

ning of the ingot R2. (400x) 
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Figure 9: Striations in middle portion of ingot R4 showed 

the length of material grown in spin-up period was different 
from that in spin-down period. (50x) 
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Figure 10: Coupled striations in grains about 6 mm from 

the tail end of ingot Rl, suggesting that both the start and 
stop of rotation cause the most significant perturbation to 
the melt flow pattern and to the heat and mass transfer. The 
average growth rate during ACRT on-time was twice as 
that during Off-time. (40x) 
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APPENDIX A 


Excerpts from PhD Thesis of Mohsen Banan 
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Chapter 4 


Results 


In this chapter the' results of compositional and microstructural analyses of the so- 
lidified ingots are presented and discussed. The results are presented in two sections: 1) 
Solidification with axial vibration, and 2) Current-induced perturbations. The Vibration 
section consists of the microstructural and compositional analyses of the ingots solidified 
with and without vibration, the vibration-induced dynamic acceleration measurements, and 
the macroscopic growth rate measurements using an interface demarcation technique. 

The Current-induced Perturbations section includes the results of solidifica- 
tion of an ingot with alternating current pulses, the current interface demarcation in an 
alloy of ln 0 . 2 Gao.sSb, and in - situ temperature measurements in the charges of GaSb and 
lno. 2 Gao sSb during passage of electric current. 

4.1 Solidification with Axial Vibration 

Several ingots with a feed composition of In 0 . 2 Ga 0 .sSb were directionally solidified 
with and without axial vibration of the ampoule. Table 3.1 shows the experimental condi- 
tions for all the runs (this table is given in the Experimental chapter). An axial temperature 
gradient of 30-35 e C/cm (measured in using a K-type thermocouple in an empty ampoule) 
and an ampoule lowering rate of 8 mm/day or 21 mm/day were used in these solidification 
experiments. The actual axial temperature in the charge would be expected to be lower 
than in an empty ampoule. Such a difference arises primarily from the higher thermal con- 


92 

PRECEDING PAGE BLANK NOT FILMED 



ductivity of a charge as compared to air. The temperature profile of the furnace used in 
our experiments was “thermally stable”, meaning that the temperature increases upward 
into the furnace. 

The freezing temperature of lno. 2 G 30 .aSb is about 690°C. Due to rejection of 
InSb, the freezing temperature decreases along the ingot. For good mixing of the melt and 
equilibrium at the freezing interface, the freezing temperature is predicted to vary from 
690°C initially to 530° near the end of growth. Comparing the furnace temperature profile 
and the range of freezing temperature, the initial solidification front was expected to be 
near the heater and adiabatic zone boundary. The freezing temperature at the final stages 
of solidification would be in the vicinity of cooler and adiabatic zone. 

4.1.1 Constitutional Supercooling in InSb-GaSb 

For solidification of Ino.jGao gSb feed composition, a groivth rate of 8 mm/day and 
an axial temperature gradient of 30-35°C/cm should avoid constitutional supercooling in 
the InSb-GaSb growth system. However, the temperature gradient in a charge is expected 
to be lower than 30-35°C/cro measured in an empty ampoule. At 21 mm/day growth rate, 
an axial temperature gradient of 30-35°C/cm or less(as in an ampoule with a charge under 
the same furnace setting) may not be large enough to avoid constitutional supercooling in 
the InSb-GaSb system. 

To show the validity of the above statement, we may use equation A. 8 describ- 
ing the conditions for avoidance of constitutional supercooling given as (see Appendix A 
for more details): 

( 4 - 1 ) 

Here m is the slope of the liquidus curve, k, is the interfacial distribution coefficient, C e 
and Cf are the total molar concentration of solute in the solid and liquid, respectively, A j 
is the mole fraction of solute in the liquid at the interface, and V is the growth rate. C c 
and Cf are calculated using the density and molecular weight data given in Appendix B 
and C. Fitting the liquidus curve for the phase diagram of InSb-GaSb we obtain: 
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(4.2) 


T = 709.5 - 93.1A'j - 96.1A', 2 

where A; is the mole fraction of InSb in the liquid. The slope of the liquidus is determined 
as: 


m = ~r = -93.1 - 192.2A, (4.3) 

a A i 

From the phase diagram, the equilibrium distribution coefficient k 0 is determined as a 
function of liquidus composition as: 

k 0 = ^ = 0.12 + 0.721 A', - 2. 37 A' 2 + 2.57A; 5 (4.4) 

A/ 

where A’ 3 is the mole fraction of solute in the solid at the equilibrium interface. Assuming 
that equilibrium prevails at the interface, the interfacial distribution coefficient is the same 
as the equilibrium distribution coefficient, i.e. k x = k 0 . 

Figure 4.1 shows plots of equation (4.1) for Gicnucai versus A / mole fraction 
of InSb in the liquid at the interface for growth rates of 8 mm /day and 21 mm/day and 
an assumed diffusion coefficient of 2xl0“ 5 cm 2 /s. The growth system is predicted to be 
stable for imposed temperature gradients above the curves and unstable for values below 
the curves for the given growth velocities. Interface breakdown is expected to appear as 
a cellular or dendritic structure resulting in axial and radial compositional fluctuations. 
The compositional variations arise from trapping of solute within the cellular structure or 
dendrites. 
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Mole Fraction InSb 


Figure 4.1: A plot of Guriticai in equation 4.1 versus mole fraction of InSb in the melt for 
InSb-GaSb system. Two growth rates were used: 8 mm/day and 21 mm/day. The system 
is predicted to be stable for an imposed temperature gradient above the curve and unstable 
for values below the curve for the given growth velocity. 
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4.1.2 Axial and Radial Compositions 

The axial and radial composition profiles of the ingots were determined using energ> 
dispersive spectrometry (EDS). The mole fraction of InSb was measured at 1 mm intervals 
along the longitudinal sections. The radial compositional profile of ingots was determined at 
1 mm intervals across the samples taken from positions along the ingot. The experimental 
EDS analysis technique is described in detail in Chapter 3. The EDS analysis calculations, 
a compositional model for good mixing in the melt, calculation of the longitudinal con- 
centration profile, and the error analysis of the EDS spectra data, are given in Appendix 
B. 

Solidification at 8 mm/day Translation Rate 

Figure 4.2 shows axial compositions versus longitudinal position of ln r Gaj_ r Sb ingots 
with feed composition of x c = 0.2. These ingots were directionally solidified at 8 mm/ day am- 
poule translation rate and 30-35°C/cm axial temperature gradient (measured in an empty 
ampoule using a K-type thermocouple). The complete mixing theory curve was calculated 
using the formulations given in Appendix B. Within experimetal error the composition pro- 
files all corresponded to good mixing in the melt for all ingots solidified, both with and 
without vibration. 

The elevated indium composition in the first to freeze portion of the ingot 
was probably due to rapid freezing following nudeation from a supercooled melt. Similar 
initial compositional variations due to delayed nudeation were also observed in directionally 
solidified Pb x Sni_ r Te [84] and In r Gai_ x Sb [6] ingots. 

Figures 4.3 and 4.4 show radial compositional profiles of ingots Dl (no vibra- 
tion) and V5 (40Hz, 0.1mm). The concentration was very uniform across the ingot. Near the 
last portion of the ingot to freeze, the compositional variation across the ingots increased. 
This radial segregation might have been due to the shape of the interface becoming more 
concave as compared to the intial sections of the ingot. Such changes in the liquid-solid 
interface shape were revealed by interface demarcation technique (details are given in the 
section on Liquid-Solid Interface Shape). A larger radial segregation is expected near the 
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Mole Fraction of InSb 



Figure 4.2: Axial composition profiles of ingots Dl(no vibration), Vl(10 Hz, 0.5 mm am- 
plitude), V2(20 Hz, 0.5 mm amplitude), and V5(40 Hz and 0.1 mm amplitude). The mole 
fraction of InSb is given versus mole fraction solidified along the ingot. All experimental 
profiles correspond to good mixing in the melt. 
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end of the ingot because the concentration change more rapidly near the end of ingot and the 
interface becomes more concave. The wavy form of the radial composition at 90% fraction 
solidified of ingot Dl (no vibration) could have been due to enhanced mixing induced by 
Maragoni convection near the top of the free melt surface. In ingot V5 (40Hz, 0.1mm), the 
radial composition variations near the end were large and asymmetrical. These variations 
could have been due to Maragoni convection at the free-melt surface, or oscillatory motion 
of the free melt surface by vibration. 



Radial Position (mm) 


Figure 4.3: Radial composition profiles of ingot Dl, which was solidified without vibration. 
Here g is the longitudinal mole fraction along the ingot. 

The axial composition profiles of ingots V3 (20Hz, 1.0 mm amplitude), V4 
(30 Hz, 0.5 mm amplitude), and V6 (100 Hz, 0.05 mm amplitude) are given in Figures 
4.5, 4.7, and 4.8. The composition profiles corresponded to good mixing in the melt. The 
experimental curve is a little higher than the theoretically calculated curve. This difference 
could be due to EDS experimental errors and to deviation of the initial feed composition 
was from 20% mole InSb. Figure 4.6 shows the radial compositional profile of ingot V3 (20 
Hz, 1.0 mm amplitude). The radial profile is uniform in composition. During solidification 
of ingots V3, V4, and V6 the furnace was shut off after partial solidification of the mell 
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in order to reveal the shape of the interface (results of interface shape are given in section 

4 . 1 . 3 ). 

The InSb-GaSb growth system is solutally stable, since the density of the 
rejected component, i.e. LnSb, is more than the preferentially incorporated component, 
GaSb. All of these ingots were directionally solidified under thermally and 6olutally stable 
conditions. Vigorous convection is not anticipated under thermally and solutally stable 
conditions. The good mixing profile in the ingots solidified with vibration was due to 
enhanced mixing by vibration of the ampoule. However, the axial composition profile of 
ingot Dl, solidified without vibration, corresponds to a profile for good mixing in the melt. 
The good mixing could have been due to convection induced by radial concentration and 
temperature gradients and/or molecular diffusion. Since the growth rate of 8 mm/day was 
very slow, there might have been sufficient time for diffusion to mix the rejected solute at 
the interface with the entire liquid (explained in the Discussion chapter). 
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Figure 4.5: Axial composition profile of ingot V3 directionally solidified at 8 mm/day 
ampoule translation rate and axial vibration of 20 Hz frequency and 1.0 mm amplitude. 
The curve corresponds to good mixing in the melt. 
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Figure 4.6: Radial composition profiles of ingot V3 which was solidified with vibration at 
20 Hz and 1.0 mm amplitude. Here g is the longitudinal mole fraction along the ingot. 
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Figure 4.7: Axial composition profiles of ingot V4 directionally solidified at 8 ram/day 
and axial vibration of 30 Hz, 0.5 mm amplitude. The curve shows the composition profile 
predicted for complete mixing of the melt. 
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Figure 4.8: Axial composition profile of ingot V6 directionally solidified at 8 mm/day and 
axial vibration of 100 Hz frequency and 0.05 mm amplitude. The line is the predicted profile 
based on complete mixing of the melt. 
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Solidification at 21 mm/day Translation Rate 


Three ingots with In 0 . 2 Gao.sSb feed composition were solidified at 21 mm/day am- 
poule translation rate and 30-35 °C/cm axial temperature gradient (measured in an empty 
ampoule). One ingot was solidified without vibration, one with axial vibration at 20 Hz 
frequency and 0.5 mm amplitude, and one with 40 Hz frequency and 0.1 mm amplitude. 

Figure 4.9 shows the axial composition profiles of ingots Fl (without vibration) 
and F2 (20 Hz, 0.5 mm) solidified at 21 mm/day. The profiles correspond to good mixing 
in the melt. Compositional fluctuations were observed after 75 to 80% mole fraction of the 
ingots had solidified. The composition fluctuations might have been due to constitutional 
supercooling and morphological breakdown. An examination of the microstructure of both 
ingots showed a change in the microstructure from multi-grain with twins to a fine grain 
structure at near where the compositional fluctuations began. After 90% mole fraction had 
solidified, the axial composition again began following a good mixing pattern. This behavior 
occurred sooner in the ingot solidified with vibration compared to the one solidified with 
vibration. Figures 4.10 and 4.11 show the radial composition profiles of ingots Fl and F2. 

Figure 4.1 (refer to the beginning of this chapter) shows the critical tempera- 
ture gradient, given in equation 4,1. in the melt plotted versus the mole fraction of InSb in 
the liquid at the interface for 8 mm/day and 21 mm/day growth rates. This figure is used 
to demonstrate the conditions for avoidance of constitutional supercooling in InSb-GaSb 
for givrn growth rates. A detail of the constitutional supercooling formulation is given in 
Appendix A and at the beginning of this chapter. An imposed axial temperature gradient 
of 30-35°C/cm (measured in an empty ampoule) was used in our solidification runs. The 
axial temperature gradient in an ampoule with a charge is expected to be lower than in an 
empty ampoule. 

For solidification runs at 8 mm/day, the temperature gradient of 30-35°C/cm 
(the actual temperature gradient measured in a charge is lower than 30-35°C/cm) is above 
the critical temperature gradient for the entire range of InSb composition in the liquid at 
the interface. This means that the interface should remain stable throughout solidification 
of the entire ingot. As mentioned in section 4.1 for ingots solidified at 8 mm/day, the 
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Figure 4.9: Axial composition profiles of In 0 . 2 Ga 0 .sSb feed composition ingots FI (no vi- 
bration) and F2 (20 Hz, 0.5 mm) directionally solidified at 21 mm/day. The continuous line 
is the theortical curve for good mixing in the melt. 
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Figure 4.10: Radial composition profiles of ingots Fl solidified at 21 mm/day translation 
rate without vibration. Note g is the mole fraction solidified along the ingot. 



Figure 4.11: Radial composition profiles of ingots F2 solidified at 21 mm/day translation 
rate with vibration of 20 Hz and 0.5 mm amplitude. Note g is the mole fraction solidified 
along the ingot. 

original page 

OF, POCm QUALITY 

105 


axial composition profiles corresponded to good mixing. No compositional fluctuations 
or other evidence of morphological breakdown were observed. In the ingots solidified at 
21 mm/day, some compositional fluctuations were observed at and beyond 75-80% mole 
fraction solidified. The compositional fluctuations might have been due to morphological 
breakdown. 

The location along the ingot that breakdown would take place could be related 
to the information given in Figure 4.1. The mole fraction InSb in the liquid at the interface 
can be determined at any given location along the ingot at the time solidification using the 
plots in Figure 4.12. Figure 4.12 shows the theoretical solute composition profile for the solid 
and liquid at the interface with equilibrium at the interface for good mixing of the melt for 
In 0 ;Gao sSb feed composition (calculation is given in Appendix B).The solute concentration 
increases in the melt ahead of interface. The mixing by free convection or \ibration-induced 
convection might not be sufficient to reduce the solute build-up. Referring to Figure 4.12. the 
liquid composition at the interface at S0% mole fraction of ingot solidified was about 0.G mole 
fraction InSb. The critical axial temperature gradient for 0.6 mole fraction of InSb in the 
liquid at 21 mm/day growth is about 100 5 C/cm. The actual growth rate with an ampoule 
lowering rate of 21 mm/day is less than 21 mm/day near the end of growth. The imposed 
axial temperature gradient during all our solidification runs was less than 30-35°C/cm . 
which is much lower than the critical temperature gradient of 100°C/cm. Consequently the 
interface should have broken down. Such morphological instability is expected to manifest 
itself as compositional fluctuations and fine grain structure, as observed in ingots Fl and 
F2. 

The interface breakdown in ingots Fl and F2 extended for 0.15 mole fraction of 
the ingot solidified and then resumed a good mixing profile again. However, the composition 
profile near the end did not follow the good mixing profile for 0.20 mole fraction InSb, since 
the solute was trapped in the broken down region. This resulted in a change of composition 
in the melt and consequently in the solid near the end. Such a recovery from morphological 
breakdown could have been due to following; after the breakdown, the solute buildup ahead 
of the interface would be reduced since the solute is trapped in the broken-down regions. 
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The rate of growth also decreases due to interaction of the growth system with the furnace 
thermal fields (the growth rate measurements are given in section 4.3). 
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Figure 4.12: Calculated axial composition profile for an ingot with a feed composition of 
lno. 2 Gao. 8 Sb with equilibrium at the interface and good mixing in the melt. The continous 
curve is the calculated InSb composition in the solid at the interface. The dashed curve is 
the InSb composition in the liquid at the interface. 
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4.1.3 Liquid-Solid Interface Shape 

Interface demarcations were generated in two ingots by lowering the heater tem- 
perature from 800°C to 770° and then abruptly translating the ampoule down by 1 mm. 
Afterward, the heater temperature was increased to 800°C again. The change of heater 
temperature would modify the thermal fields in the melt and charge, and cause fast freez- 
ing during the temperature lowering period. The ampoule was moved down by 1 mm to 
assure that the rapidly frozen demarcation would not backmelt. 

The ingots were longitudinally sectioned, mechanically polished, and chemi- 
cally etched in 1 HF: 1 HN' 03 : 1 H 2 0 for 30 sec to 2 minutes. A Longer etching time caused 
extensive cracking of the ingots. The demarcations appeared as discontinuities in the mi- 
crostructure across the ingot. The demarcated lines had a sudden change of composition. 
The interface demarcations were not clear in the first half of the ingots, even after etching 
for 1 to 2 minutes. The demarcations were more pronounced at the second half, and espe- 
cially near the end of the ingots. The less pronounced demarcations could be due to two 
reasons: the cooling period to induced fast freeze was not long enough and/or the prefer- 
ential etching of second half as compared to first half of the ingots (the end portion had 
higher concentration of InSb compared to the first half of the ingots). At the first half of 
the ingots the shape of the interface was near-concave, with a radius of curvature less than 
0.1 mm. The radius of curvature of the demarcations increased as solidification proceeded. 
In the last to freeze section of the ingots, the interface became more concave, with a radius 
of curvature of 1 mm. In these experiments no significant difference was observed in the 
demarcations’ curvature in ingots solidified with or without vibration. 

As mentioned earlier in the previous section ingots V3 (20 Hz, 1.0mm ampli- 
tude) and V6 (100 Hz, 0.05mm amplitude) were quenched after partial solidification. The 
shape of the interface in ingot V3 was convex, as shown in Figure 4.13. The first demar- 
cation was made by lowering the heater temperature similar to the procedure mentioned 
above. The second demarcation was generated by quenching the remaining melt by shutting 
off the heater. The microstructure changed across the demarcation due to the discontinuity 
in growth rate. The first demarcation was very clear. The second demarcation was not 
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pronounced. The ingot cracked in the vicinity of the second demarcation during cutting. 
The EDS measutements showed the compositional variations across the demarcations. The 
corresponding axial composition profile is given in Figure 4.14. The compositions profile of 
the directionally solidified portion corresponded to good mixing in the melt. The interface 
shape was convex with 0.25 mm radius of curvature. A convex interface is desirable for 
better grain selection, since the grains tend to grow normal to the interface. 

Figures 4.15 shows an interface demarcation generated by quenching the re- 
maining melt during solidification of ingot V6. The ingot was chemically etched using 
1 HF: 1 HN 03 : 1 H 20 for 25 sec at room temperature to reveal the microstructure and in- 
terface demarcation. The interface was concave with a radius curvature of 0.3 mm. The 
interface was wavy and asymmetric. Figure 4.16 shows the axial composition profile in ingot 
Y6. The composition suddenly changed at the demarcation location. The location of the 
interface demarcation along the ingot correspoded to freezing temperature of 600 °C. The 
temperature of 600 °C in our furnace was in the lower section of the adiabatic 2 one. It is 
expected that the interface shape would be near-concave in the boundary of adiabatic zone 
and the cooler. 
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Figure 4.13: Microstructure of ingot V3 translated at 8 mm/day with axial vibration of 
20 Hz frequency and 1.0 mm amplitude. The interface demarcations were made by rapidly 
freezing the melt. The ingot was sandblasted to reveal the microstructure and the interface 
demarcations. The interface was convex with a radius of curvature of 0.25 mm. The ingot 
diameter was 9 mm. The growth direction was from left to right. 


1 .0 


GO 


0.3 r 


£ 0.5 

o 

£3 

£ 0.4 

o 

c 0.2 h 


GQQQD 20 Hz, 1 ram amplitude. 
Complete mixing (theory)- 


1 demarcation 


l 


o 

c 

o 

3 

C 


O 


O 


O 


o 


0. 


( ppnnnnflQTt OOOO- Q^^ 


° °o o- 

J 1_ 


0.2 0.4 0.6 0.8 

Mole Fraction Solidified 


1.0 


Figure 4.14: Axial composition profile of ingot V3. 
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Figure 4.15: Microstructure of last portion of ingot V6 to freeze with a translation rate of 
8 mm/day and axial vibrations of 100 Hz frequency and 0.05 mm amplitude. The interface 
demarcations were made by rapidly freezing the melt. The ingot was chemically etched 
using 1HF:1HN0 3 :1H 2 0 for 25 sec at room temperature to reveal the microstructure and 
interface demarcation. The interface was concave with a radius of curvature of 0.3 mm. 
The diameter of the ingot was 0.9 cm. The growth direction was from left to right. 
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4.1.4 Macroscopic Growth Rate 


The macroscopic growth rate was determined for two ingots with Ino jGao gSb feed 
composition, one without vibration and with vibration at 20 Hz frequency and 0.5 mm 
amplitude. The quenching interface demarcation technique was used. Both ingots were 
directionally solidified at 9.5 mm/day ampoule lowering rate. 

In this technique the heater temperature was decreased from 800°C to 770°C 
within 10 minutes by lowering the voltage to the heater. When the heater temperature 
reached 770°C, the ampoule was abruptly moved downward by 1 mm. After the ampoule 
was moved down by 1 mm. the heater temperature was increased again to 800°C and 
maintained at that temperature. This procedure was undertaken every 24 hours during the 
entire solidification period of 15 to 16 days. After completion of a growth run. the ingot 
was removed from the ampoule. The ingot was sectioned axially and chemically polished 
using 1HF:1HN0 3 :1H 2 0 for 60 sec at room temperature, to reveal the demarcations. 

Figure 4.17 shows the measured distance between the above demarcations 
along the ingots. The measured distance versus time was fitted into a third degree polyno- 
mial. The macroscopic growth velocity along the ingot was calculated taking the derivative 
ofthe polynomial fit of the plots of distance versus time. A plot of macroscopic growth ve- 
locity versus length fraction solidified along the ingots is shown in Figure 4. IS. The growth 
velocity was initially higher than the translation rate, possibly due to the end effects as 
predicted by Sukanek [65]. As mentioned earlier in section 4.2, the axial composition profile 
of all ingots directionally solidified with and without vibration corresponded to good mix- 
ing of the melt and equilibrium at the interface. As the growth proceeds, rejection of InSb 
lowers the freezing temperature (as given by the phase diagram), resulting in variation of 
the freezing temperature along the length of the ingot. The freezing temperature was ini- 
tially about 690°C and gradually decreased to 530°C near the end of ingot. Consequently, 
solidification would never reach a steady-state condition. 
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Figure 4.17: Position of interface demarcations versus time of generation. The circles 
and triangles represent the data for ingots solidified with vibration of 10 Hz and 0.5 mm 
amplitude and without vibration, respectively. 
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Figure 4.18: Macroscopic growth velocity versus length fraction solidified. Dashed line 
represents the ampoule translation rate of 9.5 mm/day. Note that the actual macroscopic 
growth velocity was always different from the translation rate throughout the entire growth 
id. 
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4.1.5 Vibration-induced Acceleration 

Vibration of the ampoule resulted in modulation of acceleration in the direction of 
applied vibration. In our solidification experiments the growth ampoule was oscillated ax- 
ially parallel to Earth's gravity. The dynamic acceleration induced by axial oscillation of 
the ampoule was measured using an accelerometer attached to the ampoule holder. The 
accelerometer's output signal was recorded at a sampling rate of 200 Hz using a data 
acquisition system and a Zenith 280 computer. Table 4.1 shows a summary of dynamic ac- 
celerations measured during axial vibration of the ampoule at different vibration paramters. 

Figures 4.19 to 4.23 show the dynamic acceleration versus time for different 
vibrational conditions. Figure 4.19 shows acceleration fluctuations of i0.005 x 9.81 m/s 
without anv applied vibration. These fluctuations are caused by background noise in the 
laboratory. The vibration at 10 Hz frequency and 0.5 mm amplitude caused periodic accel- 
eration variations of ± 0.1 x 9.81 m/s 2 , as shown in Figure 4.20. A periodic variation of 
± 0.2 x 9.81 m/s 2 was measured for vibration at 20 Hz frequency and 0.5 mm amplitude, 
as shown in Figure 4.21. Vibration at 40 Hz frequency and 0.1 mm amplitude resulted in 
± 0.1 x 9.81 m/s 2 acceleration variations, as shown in Figure 4.22. Figure 4.23 shows the 
dvnamic acceleration fluctuations of i 0.05 x 9.81 m/s 2 measured during axial vibration 
at 100 Hz frequency and 0.05 mm amplitude. 

The acceleration data were analyzed using the Power Spectrum option of Tem- 
pleGraph software package (via Clarkson-N AS A/Lewis computational facilities). The max- 
imum g-variation was determined at its respective frequency for dynamic acceleration data 
of 20 Hz and 0.5 mm amplitude, 40 Hz and 0.1 mm amplitude, and 100 Hz and 0.05 mm 
amplitude, as shown in Figures 4.24, 4.25, and 4.26. respectively. 
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Table 4.1: Dynamic acceleration measured using an accelerometer at different vibration 
conditions. The acceleration varied periodically. Similar vibration parameters were used in 
the solidification experiments. 


Ingot 

Vibration Parameters 

Dynamic Acceleration 


(Hz). (mm) 

(g=9.Sl m/ 5 2 ) 

VI 

10.0.5 

±0.09g 

Y2 

20.0.5 

±0.1Sg 

Y3 

20.1.0 

±0.23 

Y4 

30.0.5 

±0.25 

Y5 

40.0.1 

±0.10g 

V6 

100.0.05 

±0.05g 
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Figure 4.19: Plot of dynamic acceleration versus time without application of axial vibration. 
The measured fluctuations were due to background noise in the laboratory. The acceleration 
was measured using an accelerometer. 



Figure 4.20: Plot of dynamic acceleration versus time with axial vibration at 10 Hz frequency 

and 0.5 mm amplitude. The acceleration was measured using an accelerometer. 
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Figure 4.21: Plot of dynamic acceleration versus time with axial vibration at 20 Hz frequency 
and 0.5 mm amplitude. The acceleration was measured using an accelerometer. 
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Figure 4.22: Plot of dynamic acceleration versus time with axial vibration at 40 Hz frequency 
and 0.1 mm amplitude. The acceleration was measured using an accelerometer. 
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Figure 4.23: Plots of dynamic acceleration versus time with axial vibration at 100 Hz 
frequency and 0.05 mm amplitude. The acceleration was measured using an accelerometer. 
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Figure 4.24: Power spectrum of dynamic acceleration at 20 Hz frequency and 0.5 mm 
amplitude vibration. 
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Figure 4.25: Power spectrum of dynamic acceleration at 40 Hz frequency and 0.1 mm 
amplitude vibration. 
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4.1.6 Microstructural Analysis 

In this work, the entire ingot was sectioned longitudinally as shown in Figure 4.27. 
for compositional and microstructural analyses. The compositional analysis was performed 
using a mechanically polished sample. For irucrostructural analysis, the longitudinal cross 
sections of ingots were mechanically polished and chemically etched to reveal the microstruc- 
ture (details are given in Chapter 3). 

Large number of voids (bubbles) were observed on the surface of the ingots, as 
shown in Figure 4.28. The number of bubbles was higher on the first half of the ingots as 
compared to the remaining portion. There was no correlation between the applied vibration 
and the number of bubbles on the surface of the ingots. 

The microstructure of ingots was examined using optical microscopy at 200X 
magnification and scanning electron microscopy at 200X or higher. All ingots were poly- 
crvstalline. The microstructure consisted of many straight and curved boundaries. Etching 
of curved boundaries yielded a triangular groove. The straight boundaries appeared as 
rectangular-shape grooves. Figures 4.29 and 4.30 show scanning electron microscopy mi- 
crographs of the curved and straight boundaries after chemically etching of the sampje in 
HF iHNOaihhO solution for 25 sec at room temperature. Figure 4.31 shows a scanning elec- 
tron microscopy micrographs of a straight boundary at 1000X magnification. The rough 
surfaces are possibly due to preferential etching of regions with different orientation. 

Figure 4.32 shows the microstructure of ingot Dl solidified without vibration. 
The boundaries were revealed by chemical etching the sample in 1 HF : IHNO 3 : IHjO 
solution 

for 45 sec at room temperature. A large number of small grains was observed in the first 
to freeze section of ingot Dl. This ingot consisted of a large number of curved and straight 
boundaries, with scattered small grains throughout the ingot. 

Figure 4.33 shows a photograph of the microstructure of the longitudinal sec- 
tion of ingot VI solidified with axial vibration of 10 Hz and 0.5 mm amplitude. The 
boundaries were revealed by chemically etching the sample in IHFiIHNOsTHjO solution 
for 45 sec at room temperature. Fewer small grains was observed in the first to freeze section 
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Figure 4. 27: Schematic diagram of a longitudinally sectioned ingot used for composition . 
and microstructural analyses. 
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Figure 4.28: Sandblasted In x Ga,_ r Sb ingot which was solidified with vibration at 20 He 
frequency aDd 0.5 mm amplitude. Bubbles were present on the surface of the ingot. Growth 
direction was from left to right. 
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Figure 4.29: Scanning electron micrograph of curved boundaries on a longitudinal section 
of In 0 . 2 Gao.sSb ingot V6. The boundaries were revealed by chemically etching the sample 
in 1 HF: 1 HN 03 : 1 H 20 solution for 25 sec at room temperature. The magnifications are 300x 
and 900x for the pictures on the left and on the right, respectively. The picture on the right 
is an enlargement of the area enclosed in the box in the picture on the left. The scale bar is 
10 microns and 3.3 microns for the 300X and 900X magnifications, respectively. The large 
arrow to the right indicates the growth direction. 
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Figure 4.30: Scanning electron micrograph of curved and straight boundaries on a longitu- 
dinal section of Ino ^Gao.gSb ingot V6. The boundaries were revealed by chemically etching 
in 1 HF: 1 HN 03 : 1 H 20 solution for 25 sec at room temperature. The magnifications are 200X 
and 1000X for the pictures to the left and to the right, respectively. The picture on the 
right is a 5 times enlargement of the boxed area in the picture on the left. The longest bar 
in the bottom of the picture is 100 and 20 microns for the 200X and 1000X magnifications, 
respectively. The large arrow indicates the growth direction. 
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Figure 4.31: Scanning electron micrograph of a straight boundary on a longitudinal section 
of In 0 . 2 Gao.gSb ingot V2. The surface was etched in 1 HF: 1 HN 03 : 1 H 2 0 for 25 sec at room 
temperature. The magnification is 1000X. The longest bar at the bottom of the picture is 
10 microns. The arrow indicates the growth direction. 
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Figure 4.32: Photograph of the microstructure of a longitudinal section of ingot Dl, which 
was solidified at 8 mm/day translation rate without vibration. The boundaries were revealed 
by chemically etching the sample in HF:HN 03 :H 20 solution for 45 sec at room temperature. 
The samples were cast in a resin mold for polishing. The width of the ingot is 0.9 cm. The 
growth direction was from left to right. 
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of ingot Yl as compared to ingot Dl. Interestingly, growth of large grains was initiated 
from the first to freeze section of ingot Yl. This ingot consisted of three grains initial)) 
and two grains near the end. A large number of twins initiated from the ampoule wall and 

grew inward. No small grain was observed in ingot VI. 

Figure 4.34 shows a photograph of the microstructure of a longitudinal section 
of ingot V2, which was solidified with axial vibration of 20 Hz and 0.5 mm amplitude. Ingot 
V2 consisted of many small grains in the first to freeze section. Large twinned grains were 
observed after 4 to 5 ram from the first to freeze section of ingot V2. The twins were mostly 
initiated from the ampoule wall and grew- inw-ard until reached to a curved boundary. In 
the second half of ingot Y2. there were only two grains, but heavily twinned. The twins in 
the second half were parallel to the growth direction. 

Figure 4.35 shows a photograph of the microstructure of a longitudinal section 
of ingot Y6. which was solidified at 8 mm/day with axial vibration of 100 Hz and 0.05 
mm amplitude. In the first to freeze section of ingot Y6. a large number of fine grains was 
present. The microstructure mostly consisted of grains in the first half of the ingot wit!, 
only a few twins. In the second half of the ingot, the number of grains was reduced to three 
with several twins. The last portion of this ingot was quenched to reveal the liquid-solid 
interface. The interface shape was near concave with a wavy configuration. The wavy shape 
might have been due to oscillatory perturbations caused by vibration. 
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Figure 4.33: Photograph of the microstructure of a longitudinal section of ingot VI solidified 
with axial vibration of 10 Hz and 0.5 mm amplitude. The boundaries were revealed by 
chemically etching the sample in 1HF:1HN03:1 Hj 0 solution for 45 sec at room temperature. 
The samples were cast in a resin mold for polishing. The width of the ingot is 0.9 cm. The 
growth direction was from left to right. 
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Figure 4.34: Photograph of the microstructure of a longitudinal section of ingot V2, which 
was solidified with axial vibration of 20 Hz and 0.5 mm amplitude. The structure was 
revealed by chemically etching the sample in 1 HF: 1 HN 03 : 1 H 20 for 45 sec at room temper- 
ature. The samples were cast in a resin mold for polishing. The width of the ingot is 0.9 
cm. The growth direction was from left to right. 
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Figure 4.35: Photograph of the microstructure of a longitudinal section of ingot V6, which 
was solidified with axial vibration of 100 Hz and 0.05 mm amplitude. The structure was 
revealed by chemically etching in 1HF:1HN0 3 :1H 2 0 for 45 sec at room temperature. The 
samples were cast in a resin mold for polishing. The last to freeze section shows a quenched 
liquid-solid interface. The interface was near concave. The width of the ingot is 0.9 cm. 
The growth direction was from left to right. 
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4.1.7 Statistical Analysis of Grain and Twin Boundaries 

A statistical analysis was performed on the number of curved and straight boundaries 
in different ingots. The curved and straight boundaries are reported as grain and twin 
boundaries, respectively, in this work. The number of boundaries per millimeter across an 
ingot was compared to other ingots. Also the total number of boundaries in each ingot was 
compared to the other ingots. 

The number of grain and twin boundaries was counted across the ingots along 
the longitudinally sections at 2 mm intervals. Figures 4.36 and 4.3 r show the number of 
grain and twin boundaries per millimeter across the width of the samples for the entire 
length of the ingots solidified with and without vibration. The number of grain boundaries 
was significantly lower for the ingots solidified with vibration compared to the one without 
vibration, especially in the second half of the ingots. The scatter in the data could be due to 
the presence of small grains along the examined sections. The last half of all ingots showed 
less scatter in the data as compared to the first half. Figures 4.38 shows the number of 
grain boundaries for the last half of the ingots. 

A one-sided Student's t-test with paired comparisons was used in this analysis. 
Thet-test was performed using a commercially available software package called the Number 
Cruncher Statistical System (NCSS). A detailed description of the Students's t-tesi is as 
follows: Consider two data sets which we desire to compare. Both of these sets are arranged 
in two columns: the first column lists the length fraction of the ingot solidified and the second 
column contains the straight boundaries per mm width of ingot. In order to compare them, 
the two data sets are paired based on equal length fractions of the ingot solidified. These 
pairs are listed under two columns, Cl and C2 as a separate data set. As described below 
a the two sided t-test was performed on this data set. 

1. Hypothesis proposed: The first step involved the proposal of an hypothesis, which 
was as follows: To check if the mean of the difference between Cl and C2 is greater 
than zero. 
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4.36: Number of grain boundaries per mm counted across the samples at 2 mm 
I? iIotic ingots solidified at 8 mm/day translation rate with and without vibration. 



Length (Fraction Solidified) 

Figure 4.37: Number of twin boundaries per mm counted across the samples at 2 mm 
intervals along ingots solidified at 8 mm/day translation rate with and wilhout vibration. 
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Figure 4.38: Number of grain boundaries per mm counted across the samples at 2 mm 
intervals along ingots solidified at 8 mm/day translation rate with and without vibration. 
The boundary counts are for the last half of the ingot. 
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2. Calculate the differences X,. 


JV, = Cl, - C2, 

3. Calculate the mean difference d. 

n 

where n = cumber of observations. 

4. Calculate the variance of the difference s 2 . 

a £?(*, - dl 
* ~ n - 1 


5. Calculate the standard deviation of the difference s. 

s = 


6. Calculate the standard error of the difference sj. 


C 



7. Calculate the T value t 0 . 


d 



6. Calculate the probability level. 


Probability level = P(|C| > 0) 


(4.5) 


(4.6) 


(4.7) 


( 4 > . 


( 4.9 ■ 


( 4.10 


( 4 . 11 ) 


9. We can say that Cl was greater than C2 at the probability (confidence) level obtained 
by the previous step. 

10. All the above mentioned 6teps were performed by the NCSS software. 

Table 4.2 shows the probability level that the number of grain boundaries 
(curved) per mm across the ingot, as shown in Figure 4.36, listed in the row was greater 
than the number of curved boundaries per mm in the ingots listed in the columns of the 
table. Accordingly, with 99 9c confidence, ingot Dl(no vibration) contained more grain 
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boundaries per mm than all of the ingots grown with vibration. Ingot \ 2(20Hz. 0.5mm < 
had the least number of grain boundaries. Ingot \1 (10Hz.0.5mm) and \4 (30Hz, 0.5mm) 
showed similar numbers of grain boundaries, lower than in the ingots \ 5 and \ 6. 

Table 4.3 shows the probability level that the number of twin boundaries per 
mm across the ingot, as shown in Figure 4.37, listed in the row - was greater than the number 
of twin boundaries per mm in the ingots listed in the columns. Ingot D1 (no vibration) shows 
the least number of twin boundaries per mm among all ingots. Ingot VI (10Hz,0.5mm) and 
Y5 (40Hz. 0.1mm) showed the highest number of twin boundaries among all ingots. Ingot 
V6(100Hz. 0.05mm) had the lowest number of twin boundaries among ingots solidified with 
vibration. 

Table 4.4 shows the probability level that the total number of boundaries, 
i.e. grain and twin boundaries (curved and straight), along the full length, per mm across 
the ingot listed in the row was greater than the total number of boundaries per mm in 
the ingots listed in the columns. Ingot D1 (no vibration) showed the highest number of 
total boundaries per mm among all ingots. Ingot Y2(20Hz. 0.5mm) solidified with vibration 
had the lowest number of total boundaries per mm among all ingots. The total number 
of boundaries per mm in ingot Y5(40Hz.0.1mm) was highest among ingots solidified with 
vibration and comparable with ingot Dl solidified without vibration. 

Tables 4.5. 4.6. and 4.7 show the probability level that the number of gram 
boundaries, twin boundaries, and total number of boundaries along the last half of the 
ingots listed in the row was greater than the boundaries per mm in the ingots listed in the 
columns. 

The mean of the number of grain and twin boundaries and total boundaries 
are given in Figures 4.39, 4.40 and 4.41, respectively. The error bars represent the standard 
error of the means. The standard error of the means is defined as the standard deviation 
divided by the square root of the total number of observations. The mean of the number of 
grain boundaries per mm was highest for ingot Dl (no vibration) as compared to all other 
ingots. Ingot Y2 (20Hz, 0.5mm) had the lowest mean number of grain boundaries per mm 
among all ingots. The mean number of twin boundaries per mm was lowest in ingot Dl 
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Table 4.2: Probability level that the number of grain boundaries per mm across the ingot 
listed in the row was greater than the number of grain boundaries per mm in the ingots 
listed in the columns. The grain boundaries along the full length of the ingots were compared 
in this table. These ingots were solidified at 8 mm/day translation rate with or without 
vibration of the ampoule. 


Ingot 

Dl 

VI 

V2 

V4 

V5 

V6 

Dl(no vibration) 


0.99 

0.99 

0.99 

0.99 

0.99 

Yl(10Hz.0.5mm) 



0.99 

0.55 



-1 

Y2(20Hz. 0.5mm) 







\'4(30Hz. 0.5mm) 



0.99 




Y5(40Hz. 0.1mm) 


0.99 

0.99 

0.97 



Y6(100Hz. 0.05mm) 


0.97 

0.99 

0.97 

0.59 



Table 4.3: Probability level that the number of twin boundaries per mm across the ingot 
listed in the row was greater than the number of twin boundaries per mm in the ingots lisied 
in the columns. The twin boundaries along the full length of the ingots were compared in this 
table. These ingots were solidified at 8 mm/day translation rate with or without vibration 
of the ampoule. 


Ingot 




D 

Y5 

V6 

Dl(no vibration) 

i 

H 


■ 



Vl(10Hz, 0.5mm) 

0.99 

HI 

0.99 

0.72 


0.99 

V2(20Hz, 0.5mm) 

0.62 

| 



0.75 


V4(30Hz, 0.5mm) 

0.95 


0.95 



0.93 

V5(40Hz, 0.1mm) 

0.99 

0.62 

0.99 

0.59 


0.99 

Y6(100Hz. 0.05mm) 

0.52 







137 

























Table 4.4: Probability level that the total number of boundaries per mm across the full 
length of the ingot listed in the row was greater than the total number of boundaries per 
mm in the ingots listed in the columns. The total number of boundaries along the full 
length of the ingots were compared in this table. These ingots were solidified at 8 mm/da\ 
translation rate with or without vibration of the ampoule. 


Ingot 

Dl 

VI 

V2 

V4 

V5 

V6 

Dl(no vibration) 


0.93 

0.99 

0.99 

0.68 

0.99 

Vl( 10Hz.0.5mm) 



0.99 

0.71 


0.99 

Y2(20Hz. 0.5mm) 







Y4(30Hz.0.5mm) 



0.99 



0.99 

Y5(40Hz. 0.1mm) 


0.87 

0.99 

0.76 


0.99 

Y6(100Hz. 0.05mm) 



0.99 





Table 4.5: Probability level that the number of grain boundaries per mm across the ingot 
listed in the row was greater than the number of grain boundaries per mm in the ingot? 
listed in the columns. The grain boundaries along the last half of the ingots, as shown 
in Figure 4. 38. were compared in this table. These ingots were solidified at 8 mm/day 
translation rate with or without vibration of the ampoule. 


Ingot 

Dl 

VI 

V2 

V4 

V5 

V6 

Dl(no vibration) 


0.99 

0.99 

0.99 

0.99 

0.99 

Vl(10Hz, 0.5mm) 



0.99 

0.55 



V2(20Hz, 0.5mm) 


■ 


BH 



V4(30Hz, 0.5mm) 



0.99 


Hi 


V5(40Hz. 0.1mm) 



mm 

UU 

■ 

0.99 

Y6(100Hz, 0.05mm) 


g§ 
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Table 4.6: Probability level that the number of twin boundaries per mm across the ingot 
listed in the row was greater than the number of twin boundaries per mm in the ingots 
listed in the columns. The twin boundaries along the last half of the ingots were compared 
in this table. These ingots were solidified at 8 mm/dav translation rate with or without 
vibration of the ampoule. 


Ingot 

D1 

VI 

V2 

V4 

V5 

V6 

Dl(no vibration) 







Vl( 10Hz, 0.5mm) 

0.99 


0.96 

0.73 


0.99 

Y2(20Hz. 0.5mm) 

0.92 




0.58 

0.54 

\'4(30Hz. 0.5mm) 

0.99 





0.86 

V5(40Hz. 0.1mm) 

0.99 

0.98 

0.99 

0.60 


0.99 

Y6( 100Hz. 0.05mm) 

0.88 







Table 4.7: Probability level that the total number of boundaries per mm across the ingot 
listed in the row was greater than the total number of boundaries per mm in the ingots 
listed in the columns. The total number of boundaries along the last half of the ingots were 
compared in this table. These ingots were solidified at 8 mm/day translation rate with or 
without vibration of the ampoule. 


Ingot 

Dl 

VI 

V2 

V4 

V5 

V6 

Dl(no vibration) 



0.99 

0.98 


0.96 

Vl(10Hz, 0.5mm) 

0.92 


0.99 

0.71 


0.98 

V2(20Hz, 0.5mm) 







V4(30Hz, 0.5mm) 



0.95 



0.89 

V5(40Hz, 0.1mm) 


0.99 

0.99 

0.72 


0.99 

V6(100Hz, 0.05mm) 



0.99 
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(no vibration) among all ingots, especially in the last half of the ingots, as shown in Figure 
4.43. The total number of boundaries was lowest in ingot V2(20Hz.0.5ram) and highest in 
ingot V5(40Hz, 0.1mm). The high number of total boundaries in ingot V5 was due to large 
number of twin boundaries. 

E 
E 

\ 
i r. 

■Z 1-5 

ec 
"C 

c 

o 
c 

B 1.0 

f: 
u 

C 

u 

c 

= 0.5 

c: 
c 
2 

0.0 

Figure 4.39: Mean number of gTain boundaries per mm width counted across the samples 
at 2 mm intervals along full length of the ingots solidified with and without \ibration a, 8 
mm/dav translation rate. The error bars represent the standard error of the mean. 

Figures 4.44 and 4.45 show the mean number of grain and twin boundaries per 
mm of the ingots solidified with and without vibration versus the vibration-induced dynamic 
acceleration. The acceleration was measured using an accelerometer during application of 
vibration. The ± accelerations represent the maximum and minimum acceleration values 
measured at specific vibration conditions. The dynamic acceleration variations followed a 
sinusoidal harmonic pattern (refer to Acceleration Measurement section in the Results chap- 
ter for more detail). The mean number of grain boundaries was highest for lg(gravitational 
acceleration without vibration). The mean number of grain boundaries decreased as the 
dynamic acceleration increased up to ±0.18, for ingot V2(20Hz, 0.5mm). A higher dynamic 
acceleration of ±0.25 (ingot V4 - 30 Hz, 0.5mm) showed an increase in the number of grains 
compared to the other ingots grown with vibration; still the number of grain boundaries 



D1 VI V2 V4 V5 V6 
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No vibration 
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20Hz. 0.5mm 
30Hz. 0.5mm 
40Hz. 0.1mm 
100Hz. 0.05mm 


*5 1.2 
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Figure 4.40: Mean number of twin boundaries per mm width counted across the samples 
at 2 mm intervals along full length of the ingots solidified with and without vibration a. 8 
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Figure 4.41: Mean number of the total boundaries per mm width counted across the samples 
at 2 mm intervals along the full length of the ingots solidified with and without vibration 
at 8 mm/dav translation rate. The error bars represent the standard error of the mean. 



Ingots 


Figure 4.42: Mean number of grain boundaries per mm width counted across the samples 
at 2 mm intervals along last half of the ingots solidified with and without vibration at £ 
mm/day translation rate. The error bars represent the standard error of the mean. 



Ingots 


Figure 4.43: Mean number of twin boundaries per mm width counted across the samples 
at 2 mm intervals along last half of the ingots solidified with and without vibration at 8 
mm/day translation rate. The error bars represent the standard error of the mean. 
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was lower than in the ingot grown without vibration. 

The mean number of twin boundaries in the ingots did not follow a well- 
defined pattern with respect to the dynamic acceleration. However, the mean number of 
twin boundaries was lowest for ingot Dl(no vibration) as compared to ingots grown with 
vibration. 
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Dynamic Acceleration x 9.81 m/s 2 


Figure 4.44: Mean number of grain boundaries per mm of the ingots solidified with and 
without vibration at 8 mm/day translation rate. The grain boundary counts are plotted 
versus the dynamic acceleration induced by vibration during growth. For example, the 
2 ero acceleration is for run Dl without vibration and ±0.05g represents the minimum and 
maximum dynamic accelerations during growth of ingot V6. 
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Figure 4.45: Mean number of twin boundaries per mm of the ingots solidified with and 
without vibration at 8 mm/day translation rate. The twin boundary counts are plotted 
versus the dynamic acceleration induced by vibration during different growth runs. For 
example, the zero acceleration is for run Dl without vibration and the ±0.05g represents 
the minimum and maximum dynamic acceleration during growth of ingot V6. 
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4.2 Current-induced Perturbations Studies 


In this section the results of a solidification run with passage of alternating current 
through the ampoule, current interface demarcation of Te-doped Ino.iGao.gSb ingot, and 
in-situ temperature measurements in the melt and in the vicinity of liquid-solid interface of 
GaSb and InGaSb charge are presented and discussed. 

4.2.1 Solidification with Alternating Current 

An ingot with ln 0 . 2 Gao.sSb initial feed composition was directionally solidified at 8 
mm/dav ampoule translation rate and 25-30 c ’C/cm axial temperature gradient (measured in 
an empty ampoule) with application of alternating 15 amp/cm 2 current pulses. The passage 
of current through the growth system was started after 72 hours of ampoule lowering. The 
current was applied alternatively with 25 sec ( 4 ) and 25 sec (-) polarities during growth. 

The 25 sec puse duration was similar to one of the ACRT experiments performed by Gra\ 

[20' at Clarkson. The original objective of this experiment was to compare an ingot solidified 
with ACRT at 25 sec cycle time to an ingot grown with current pulses of 25 sec duration. 

The common point between the ACRT and current-induced growth was the possibility of 
backmelting and regrowth behavior. The magnitude of periodic growth for InSb-GaSb 
under application of ACRT or electric current is not available. 

Figure 4.46 shows a photograph of the longitudinal section of this ingot. The 
sample was sandblasted to reveal the microstructure. This ingot had a much finer grain 
structure with microcracks than the ingot shown in Figure 4.47, solidified without current 
pulses under otherwise identical growth conditions. 

The axial and radial compositional profiles of ingots were measured using EDS 
(details described in Experimental section). Axial composition profiles of the above two in- 
gots are shown in Figure 4.48 and 4.49. The cross-sectional composition profiles of these 
ingots are shown in Figures 4.50 and 4.51. The ingot solidified without current shows an 
axial composition profile corresponding to good mixing in the melt and uniform radial com- 
position at different positions in the ingot. On the other hand, the ingot Cl solidified with 
current shows significant axial and radial compositional variations. The radial composition 
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in the initial section, up to about g=0.12, of ingot Cl is fairly uniform. This section of the 
ingot was directionally solidified without current. 

The microcracks in ingot Cl were possibly due to stress induced by composi- 
tional variations. The fine grain structure, microcracks, and compositional fluctuations in 
ingot Cl, which was solidified with current pulses, are indicative of morphological breakdown 
caused by passage of current. The criterion for occurrence of constitutional supercooling in 
a binary alloy with and without application of electric current is formulated and given in 


Appendix A. 




Figure 4.46: Longitudinal section of an ingot with feed composition In 0 . 2 Ga 0 .sSb frozen at 
a translation rate of 8 mm/day with application of alternating 15 amp/cm^ current. The 
ingot was sandblasted to reveal the microstructure. The microstructure exhibits a fine grain 
structure with microcracks. The growth direction was from left to right. 

From Appendix A, a heat balance at the interface can be written with and 
without current applied: 

(4.12) 

(4.13) 
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k t G a - k t Gi = VncH (No Current ) 
k a G t - k t Gi = (V^ + Vc)U + *■/ (With Current ) 
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Figure 4.47: Longitudinal section of an ingot with feed composition Ino.jGao.sSb direction- 
ally frozen at an ampoule translation rate of 8 mm/day without application of current. The 
ingot was chemically etched in 1HF:1HN03:1H20 solution for 20 sec at room temperature 
to reveal the microstructure. This ingot consists of several twinned grains. The growth 
direction was from left to right. 
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Figure 4.4S: Axial compositioDal profile. 1 ~ T an ingot directionally solidified without app.i- 
cation of current. The profile corresponds ?o good mixing in the melt. 



Figure 4.49. Axial compositional profile of an ingot solidified with alternating 15 amp/cm 2 
current pulses. Current was applied after initial 4 days of ampoule lowering corresponding. 
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Figure 4.50: Radial compositional profiles of the ingot directionally solidified without ap- 
plication of current. Here g is the longitudinal mole fraction down the ingot. 



Figure 4.51: Radial compositional profiles of the ingot Cl solidified with application oj 
current. Here g is the axial position expressed as mole fraction of the ingot solidified. 
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where V c is the interface velocity without current applied, G t and Gi are the axial tem- 
perature gradients in the solid and liquid at the interface, respecti\el\. x is the Peltier 
coefficient, R is the heat of fusion per unit volume, and I is the current density. The net 
velocity during passage of current is (V^ + V c ). The perturbation in temperature gradients 
is negligible at the instant of application of current. At that instant the current-induced 
change in growth rate is, from the above equations: 


\ 


r 

c — 


-7T J 

TT 


(4.14) 


Equations describing the conditions for avoidance of constitutional supercooling are gi\en 
as (see Appendix A for more details): 


A / m G c . _ . . - , 

(.Vo Current) Gi cr ,t,ca! > — nc(*t ~ 1J 


(4.15 


At the first instant when current is turned on: 


Guinea! > Y ){k, ~ 1)" l (pI} ' i 


(4. 1C' 


where tu is the slope of liquidus curve. k t is the interfacial distribution coefficient determined 
from phase diagram, C c and C j are the total molar concentration of solute in the solid and ir. 
the liquid at the interface, respectively. V is the electromobility, p is the electrical resistivity, 
and 1 is the current density. 

The above equations were used to estimate the critical temperature gradient 

Gicnticai as a function of composition for In*Gaj_ r Sb. 

In this work Ino. 2 Gao. 8 Sb feed composition was used for all directional solidi- 
fication experiments both with and without current pulsations. Fitting the liquidus curve 
for the phase diagram of InSb-GaSb, shown in Figure 4.52, we obtain: 


T = 709.5 - 93. 1A',-96.1A? (4.17) 

where A'/ is the mole fraction of InSb. The slope of the liquidus is determined by differen- 
tiation of the above equation as: 
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Figure 4.53: Equilibrium distribution coefficient k c versus mole fraction A'/ liquidus con 
centration of InSb for InSb-GaSb system. The circles are the data taken from the phase 
diagram. The line is a polynomial fit given as k 0 = = 0.12 + 0. 721A';-2.37A; J + 2. 5. A/ 3 . 


152 


ORIGINAL PAGE JR 
OF POOR (Jual/A 





m = —r - -93.1 — 192.2 A, 
cA i 


(4.18) 


From the phase diagram, the equilibrium distribution coefficient k 0 is determined as a 
function of the liquidus composition of InSb, as shown in Figure 4.52: 


k e = = 0.12 + 0.721 A'< - 2.37A? + 2.57A7* (4.19) 

A/ 

where A', is the mole fraction of solute in the solid at the interface. Assuming that equi- 
librium prevails at the interface, the interfacia] distribution coefficient is the same as the 
equilibrium distribution coefficient, i.e., k , = k 0 . 

The thermophysical properties of IniGaj_ r Sb are not known versus compo- 
sition. The properties of InSb and GaSb are fairly well determined and are given in Ap- 
pendix C. Here in our calculations we used properties for GaSb. such as Peltier coefficient 
( ~ CzS: = 0 -06 volt), latent heat of fusion (Hc a s!. = 1300 J/cm 3 ). and resistivity (p=0.00, 
0. cn - . For estimation of the critical G; the growth rate without application of current was 
set to the ampoule lowering rate of 8 mm/day. which is 9.25 x 10 -6 cm/s. 

In principle, electromigration in a ternary system, such as molten In-Ga-Sb. 
cannot be described correctly using a formulation for a binary system. In molten In-Ga- 
Sb. the stochiometric {In + Ga)/Sb will not be maintained in the presence of diffusion 
and electromigration. However, there is no information available on electromigration or 
diffusion in molten ln-Ga-Sb in the literaiure. In the absence of data on ternary transport 
properties, we assumed that the stoichiometry is maintained as a pseudo- binary mixture 
of InSb and GaSb, recognizing that the results are only approximate. InSb is considered 
to be the solute. The electromobility coefficient U is assumed to be the relative migration 
between InSb and GaSb in the presence of an applied electric field. Values of D = 2xl0 _!> 
cm 2 /s and l’ = lxl0 -4 cm 2 /s.volt were used. These are within the range of the values for 
metallic melts [80]. 

In Figure 4.54 the critical temperature gradient Gi is plolted versus the mole 
fraction of InSb for the InSb-GaSb system. As noted above, this plot is not valid to quan- 
tify the effect of electromigration in the InSb-GaSb system. It is used only to qualitatively 
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demonstrate the effect of electromigration in a pseudo-binary system. Note that the electro- 
migration induced by a current density of 15 amp/cm 2 is not predicted to show a significant 
effect on the critical temperature gradient, especially at lower values of mole fraction. 

Figure 4.55 show-s a plot of critical temperature gradient including the Peltier 
effect, induced by application of *15 amp/c: * current density. The electromigration term 
was neglected. The negative polarity current results in Peltier cooling at the interface, 
therefore momentarily increasing the rate of growth. Due to application of *15 amp/cm 2 
current, the instantaneous current-induced growth increase is estimated to be: 

Y c = — — = 0.0007 cm/s (4.20 

H 

Comparing the growth rates with and without current-induced Peltier cooling we find: 

0 °°^2_ = to (4.21. 

\ r: 0.000009 

Thus at the initial instant of Peltier cooling, the interfacial velocity is estimated to be 
about 70 times higher than before the current was applied. During application of current 
the thermal field in the melt and solid changes. Consequently the interfacial velocity shoun. 
decrease from its initial sudden rise and approach the rate before application of current. 

Brush et al. [91] calculated the interfacial velocity for solidification of InSb 
during application of electric current (discussed in detail in the section on Current-induced 
Perturbations). For repeated current pulses of 9.5 amp/cm 5 on for 20 sec and off for 40 sec. 
as shown in Figure 4.56, the interfacial velocity was predicted to suddenly increase from 
10 micron/sec ampoule lowering rate to 18.5 micron/sec. Due to thermal relaxation, the 
interface velocity was predicted gradually to decrease and then suddenly to fall to below 
the ampoule translation rate when the current is turned off. 

The estimated 70 fold increase in freezing rate caused by application of current 
during solidification of InSb-GaSb alloy system would require a very high temperature 
gradient to avoid constitutional supercooling. In our work an imposed temperature gradient 
of less than 25-30 °C/cm was used in the directional solidification of InSb-GaSb ingots 
(discussed in detail in the Experimental section). Although this is below the gradient 
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Figure 4.54: Plot of critical temperature gradient versus liquid composition for JnSb-GaSb 
system with current-induced electromigration. These parameters were used; D = 2xl0" 5 
cm 2 /s, H C oSt = 1300 J/cm 3 , p= 0.001 tt. cm, U=lxl0 -4 cm 2 /s.volt and V nc = 9.25 x 10~ 6 
cm/s (8 mm/day). Only the electromigartion term was considered. The Peltier term was 
neglected. The curves represent the cases for 1) 1=0 no current-induced effect, 2) I= + 15 
amp/cm 2 current density, and 3) I=-15 amp/cm 2 current density. Above the curves it is 
6table condition and below the curves, it is unstable. 
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Figure 4.55: A semi-log plot of critical temperature gradient versus liquid composition 
for InSb-GaSb system with current-induced Peltier cooling. These parameters were used 
D=2xl0 _s cm 2 /s, Hc a st = 1300 J/cm 3 , p=0.001 fi.cm, \ , nc =9.25xl0- 6 cm/s (8 mm/day). 
and jr Ca5fc =0.06 volt. Only the Peltier term was considered. The electromigration term was 
neglected. The curves represent the cases for; 1) 1=0 no current-induced effect, 2) I = -15 
arop/cm 2 current density (Peltier cooling effect). Above the curves, interface is stable and 
below the curve, it is unstable. 
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Figure 4.56: Plot of calculated interfacial velocity as a function of time for LnSb during 
periodic application of 9.5 imp/cm 5 current for 20 sec on and 40 6ec off [91]. 


157 


ORIGIN# w v ^GE IS 
OF POOP QUALITY 


required to avoid constitutional supercooling, the interface may not break down during a 
pulse. Instead new grains and twins may nucleate in the momentarily supercooled melt, 
before the interface can change shape. This may result in formation of fine grain structure 
and cellular growth. In our growth experiment with applied current the microstructure of 
the ingot consisted of fine grains with significant compositional fluctuations. 
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4.2.2 Current Interface Demarcation of Te-doped InojGao gSb 

A current interface demarcation experiment was performed with Inj.Gaj_j.Sb doped 
with tellurium. A 7 cm long by 9 mm diameter charge of Ino.zGao.gSb was allojed in the 
rocking furnace for 9 hours at 820°C. The charge was placed in an electroded ampoule and 
1000 ppm by weight of six-nines purity tellurium was added to the charge. The ampouie 
was purged with argon and sealed under a vacuum of 10 -6 torr. The ampoule was placed in 
the Bridgman-Stockbarger furnace with setting? of 800°C and 475°C for heater and cooler, 
respectively, and a 5 cm long adiabatic zone. The charge was allowed to melt down aand 
left over night. The ampoule was occasionally shaken manual!) for homogenization Oi Te 
dopant in the melt. The ampoule was lowered at 8 mm/day. 

After solidification, the resulting ingot was sectioned longitudinally, cast in 
a resin mold, and mechanically polished. These samples were chemically etched in 1HF . 

1HN0 3 : lKMn0 4 solution for 60-80 sec at room temperature, and rinsed in de-ionized 
waier and methanol. The samples were examined using darkfield and Nomarski optical ms- 
croscopv and scanning electron microscopy. Current-induced growth rate \ariations man- 
ifested themselves in the ingot as Te concentration bands known as pulse striations [2V. 

The Te-rich regions were preferentially etched, resulting in bands with different topograph) 
compared to the regions solidified without current. 

A series of 10 amp current pulses with 10 sec on-time and 300 sec off-time 
duration were tried. The current was passed from solid ( + ) to melt (-). The current 
pulses were passed during the last 5 days of growth. The ingot had a poor microstructure 
with microcracks and many small grains. Striations were oberved in some grains and not 
in others. Consequently, it was difficult to trace these striations to determine the overall 
shape of the liquid-solid interface at the time of a pulsation. 

Figure 4.57 shows a scanning electron micrograph of pulsed striations in Ino 2 Ga 0 sSb 
feed composition ingot doped with tellurium. These striations were generated by passage of 
10 amp current (15.7 amp/cm 2 ) for 10 sec followed by no current for 300 sec. The current 
was passed from solid ( + ) to melt (-). The groved boundary is a grain boundary with 
the interface demarcations crossing the boundary. Interestingly, the demarcations were ir- 

CRK3IN I CAGE IS 
OF POOR QUALITY 


159 



regular in the vicinity of some boundaries. Figures 4.58 and 4.59 show photomicrographs 
of striations across straight and curved boundaries. These photomicrographs were taken 
using Nomarski optical microscopy with polarized light. A closer examination of the curved 
boundary in Figure 4.58 shows that the boundary was wigglv. The frequency of striations 
corresponded to the frequency of the wiggle. Some irregularities in demarcated regions were 
observed in the vicinity of the wiggled, curved boundary. 

Figure 4.60 shows a photomicrograph of striations crossing a twin boundary 
in Te-doped lno. 2 Gao.gSb generated by passage of 10 amp current from solid ( + ) to melt (-) 
for 10 sec followed by no current for 300 sec. The striations were revealed by etching of the 
sample in lHF:lHN0 3 .lKMn0 3 solution for 1 min at room temperature. The photomicro- 
graph was taken using Nomarski optical microscopy with polarized light. The demarcations 
were regular across the twin boundaries. No sign of instability or irregularity was obserceu 
in the demarcation across the twin boundaries. A slightly shift of demarcation position was 
observed with twinned regions, as seen in Figure 4.60. 
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Figure 4.57: Scanning electron microgrpah of striations in a Te-doped In 0 . 2 Gac.sSb ingot. 
A 10 amp current was passed from solid ( + ) to melt (-) for 10 sec followed by 300 sec 
without current. The sample was etched in lHF:lHN 03 :lKMn 03 solution for 1 min at the 
room temperature. The magnification is 500X. The long scale bar, pointed at with a small 
arrow, is 10 micron. Direction of growth was as shown by the long arrow in the right upper 
right-hand corner of the picture. 
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Figure 4.58: PhotomicrogTpah of striations in a Te-doped Ino.jGao.gSb ingot generated by 
passage of 10 amp current from solid ( + ) to melt (*) for 10 sec followed by no current for 300 
sec. The striations were revealed by etching the sample in lHF:lHN 03 -lKMn 03 solution 
for 1 min at room temperature. This photomicrograph was taken using Nomarski optical 
microscopy with polarized light. Note that the striations cross one curved and one straight 
boundary. The magnification was 100X. The direction of growth was from bottom to top 
of the picture. 
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Figure 4.59: Photomicrogrpah of irregular striations in a Te-doped In 0 . 2 Ga 0 .gSb ingot gener- 
ated by passage of 10 amp current from solid (+) to melt (-) for 10 sec followed by no current 
for 300 sec. The striations were revealed by etching the sample in lHF:lHN 03 .lKMn 03 
solution for 1 min at room temperature. This photomicrograph was taken using Nomarski 
optical microscopy with polarized light. Note the irregularity in the striations near the 
boundary. The magnification is 200X. The direction of growth was bottom to top of the 
picture. 
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Figure 4.60: Photomicrograph of striations crossing a twin boundary in Te-doped 

In 0 2Ga 0 .gSb generated by passage of 10 amp current from solid ( + ) to melt (-) for 10 
sec followed by no current for 300 sec. The striations were revealed by etching of the sam- 
ple in lHF:lHNC>3.lKMn03 solution for 1 min at room temperature. The Photomicrograph 
was taken using Nomarski optical microscopy with polarized light. The direction of growth 
was from bottom to top of the picture, magnification was 225X. 
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4.2.3 In-silu Temperature Measurements with Current 

Passage of current through a melt-solid system results in Peltier heating or cooling 
at the interface, Joule beating, and Thomson heating or cooling in the solid and meh. 
These thermoelectric effects cause transient thermal perturbations in the melt and solid 
The magnitude of the thermal perturbations depends on the magnitude and direction of the 
current(for reversible Peltier and Thomson heating or cooling). In semiconductor materials, 
the Peltier effect depends on the electronic behavior of materials. For example, in n-type 
materials such as InSb (intrinsicly n-type), passage of current from melt( + ) to solid(-) 
results in Peltier heating. By changing the direction of current. Peltier cooling results. 
Whereas in p-tvpe materials such as GaSb. the effect of application of current at the me.t 
- solid interface is opposite to n-type materials. 

An alloy of InSb-GaSb changes from p-type to n-type at about 4591 mole 
fraction InSb concentration in the solid [36'. Generating current-induced Peltier heating 
or cooling at the melt-solid interface in the InSb-GaSb system depends on the electronic 
behavior of the system. Due to rejection of InSb at the interface, the growing ingot contain 
less InSb initially and more near the end. It means that the materials changes from p- 
type to n-type somewhere along the ingot. If continuous Peltier cooling pulses are needed 
to demarcate the interface, the direction of current must be switched from one polarity 
to another when the ingot changes from p-type to n-type. Te-doping also changes the 
electronic behavior of the InSb-GaSb system. 

In this work, we tried to determine the dependence of Peltier effect on the 
current polarity in GaSb system. Also we tried to measure the thermal perturbations 
induced in the melt and in the vicinity of the liquid-solid interface in GaSb. These results 
were used to separate the effects of Joule heating, Thomson heating and cooling, and Peltier 
heating and cooling in GaSb. 

In this section we report results of in-situ temperature measurements in melts 
of In r Ga,_ r Sb and in GaSb during passage of electric current pulses. The temperature 
measurement technique is given in Chapter 3. In brief, an ungrounded K-type sheathed 
thermocouple was placed in a charge of in a vacuum sealed ampoule. The charge was 
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contacted on top and bottom by two molybdenum-graphite electrodes used for passage of 
electrical current. The available programmable current generator could supply up to 10 amp 
current. Current amplitudes of 2 to 10 amp were tried. Here we report the temperature 
measurements for 10 amp current. The thermocouple was connected to a thermocouple 
board of a data acquisition system and a Zenith 248 computer. Pulse durations of 5. 30. 

and 60 sec were used in these experiments. 

Overall, seven individual ampoules were prepared for temperature measure- 
ments. Only three experiments were successful. The other four were terminated due to 
failure of the thermocouple sheath or by leakage of melt from the bottom electrode, result- 
ing in loss of electrical contact between the electrodes and charge. 

Temperature Measurements in Molten Inc.:GaogSb 

Several temperature measurements were performed in a molten charge of prealloyed 
In 0 ;Ga 0 $Sb. The charge was 7 cm long and situated in a quartz ampoule of 9 mm inner 
diameter and 11 mm outer diameter. An ungrounded K-type thermocouple with 0.41 mm 
diameter 304 stainless steel sheath was placed 3.5 cm into the charge from bottom. In 
these measurements a strip-chart recorder and digital thermometer were used to record 
and monitor the thermocouple readings in the charge. The entire charge was situated in 
the heater section of the furnace. The healer and cooler settings were 820^0 and 450 : C. 
respectively. After melting the charge, the thermocouple reading in the charge was allowed 
to reach steady-state before application of current. 

Figure 4.61 shows the thermocouple readings in the melt of In 0 . 2 Gao.gSb dur- 
ing periodic passage of 10 amp current for 60 sec followed by no current for 30 off. The 
thermocouple reading was initially at 815. 7°C before application of current. The tempera- 
ture periodically varied between 815. 7°C to 817.2°C. The temperature in the melt increased 
gradually during current passage and decayed when it was off. 

Figure 4.62 shows thermocouple readings in the In 0 . 2 Gao.gSb melt during peri- 
odic passage of 10 amp alternating current pulses for 30 s with positive polarity and 30 with 
negative polarity. The thermocouple reading was initially at 808. 9°C before application of 
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current. The temperature of the melt increased due to a new steady state of about 81 1.5 C. 
The temperature remained fairly constant during alternating current pulses. Due to the 
low resolution of the strip-chart recorder, the periodic Thomson cooling and heating were 
not detectable. 
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Figure 4.61: Thermocouple readings in molten Ino. 2 Ga 0 .sSb during passage of 10 amp cur- 
rent for 60 sec followed by no current for 30 sec. The current was passed from the top 
electrode to the bottom electrode. The tip of the thermocouple was situated 3.5 cm into 
the 7 cm long charge. 
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Temperature Measurements in GaSb 

Several in-sift; temperature measurements were performed in GaSb during passage 
of current. The objective of these experiments was to measure the current-induced thermal 
perturbations in GaSb due to Joule heating, Thomson heating or cooling and Peltier effect. 
GaSb was used instead of Ino. 2 Gao. 8 Sb. The freezing temperature of InSb-GaSb vanes 
along the ingot. It was difficult to determine the exact position of the thermocouple tip 
with respect to the freezing temperature of the InSb-GaSb charge. By knowing the freezing 
temperature of GaSb, it was possible to position the tip of thermocouple junction in the 
vicinity of the interface. These measurements were performed using a K-type thermocouple 
with stainless steel sheath placed in a charge of GaSb. The tip of the thermocouple wa? 
positioned 3.5 cm into a 7 cm long GaSb charge. The success of these experiments was 
limited by the failure of the thermocouple in the melt. The melt soon damaged the sheath 
and the thermocouple junction. These measurements were performed with heater and cooler 
settings of 800°C and 475°C. respectively, a 5 cm insulation layer, and a stationary ampoule. 

Figure 4.63 shows the thermocouple readings in the GaSb melt during passage 
of 10 amp current alternating for 30 sec (4-) polarity and 30 sec (-) polarity In recording T 1 . 
the molten charge was situated in the middle of the heater. Recording T2 was made aFer 
the ampoule containing the charge and thermocouple was moved 1.5 cm down from where 
recording Tl had been taken. The axial temperature gradient was higher at the position 
where T2 was taken than where Tl was taken. In Tl. the passage of current was initiated 
after 20 sec. The thermocouple reading showed a temperature rise from 780. 5°C to 782.2°C . 
and then the thermocouple reading remained at about 782°C throughout the remaining of 
pulsations, indicating a new near-equilibrium was established. Small temperature fluctua- 
tions were observed during the change of direction of current due to the contribution of the 
Thomson effect. 

In T2 the temperature increased from 774°C to 776.8° and afterward varied 
periodically by 0.7°C during positive and negative alternating current. 

In both Tl and T2, the temperature rise was due to Thomson heating and 
Joule heating. The periodic temperature rise and fall were due to the contribution of Thom- 
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son heating and cooling as the direction of current was changed. The periodic Thomson 
heating and cooling was less pronounced in Tl as compared to T2. Such a difference was 
probably due to the higher temperature gradient in the molten charge where the profile T2 
was recorded as compared to the Tl. The presence of liquid -solid interface was unlikeh in 
these measurements. If any solid phase was present, the liquid-solid interface would ha\e 
been about 3 cm away from thermocouple tip in recording T2. Peltier cooling possibh had 
some effect on these thermal perturbations. 

The Thomson effect is a reversible phenomenon depending on the direction of 

current. The Thomson heat is expressed as: 

where the temperature gradient. I is the current density . and r is the Thomson 

coefficient. Joule heating is expressed as: 

Q : = pl‘ H23; 

where p is the electrical resistivity and / is the current density. Joule heating is independent 
of the direction of applied current. 

A single pulse measurement was performed at the same position where profile 
T2 was taken. Figure 4.64 shows the thermocouple reading in the GaSb melt during passage 
of 10 amp current. The pulse was passed for 60 sec from the top electrode (-) to bottom 
electrode ( + ), and then vice versa. The thermocouple readings when the pulse was passed 
from top electrode (-) to bottom electrode ( + ) showed a 0.6°C temperature difference 
due to Thomson effect by passing different polarity current. In these measurements, the 
temperature rise in the melt was due to Joule heating and Thomson heating when the 
current was passed from lop(-) to bottom( + ). The reversed polarity resulted in Joule 
heating and Thomson cooling in the melt. In these measurements, the presence of liquid- 
solid interface was unlikely. 

Table 4.8 shows the power per unit volume generated in a GaSb melt for 
15.7 amp/cnr 2 current density due to Joule heating and Thomson effect. Two temperature 
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Figure 4.63: Thermocouple reading? in the GaSb melt during passage of 10 amp currer,. 
alternating for 30 sec ( + ) polarity and 30 sec (-) polarity. The Passage of current was 
initiated after 20 sec on the time scale. Recording Tl and T2 were taken at position? 
with low and high axial temperature gradients, respectively. The temperature increase was 
due to Joule heating and Thomson effects. Larger periodic thermal perturbations were 
measured T2 as compared to Tl due to the periodic contribution of Thomson heating and 
cooling because of larger temperature gradient in T2. 
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Figure 4.64: Thermocouple readings in a GaSb melt during passage of 10 amp current 
for 60 sec passed from top electrode(-) to bottom electrode(+), and the vice versa. The 
thermocouple tip was in the same position as T2 in Figure 4.63. The temperature increase 
was due to Joule beating and Thomson effects. The temperature difference of 0.6 °C was 
due to contribution of Thomson heating or cooling. 
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gradients were used in these calculations; 2°C/cm for low - temperature gradient region 
in the melt, such that in the middle of the beater, and 20°C/cm a for high temperature 
gradient region, as in the vicinity of the liquid-solid interface. For 10 amp current, the power 
generated by Joule heat is about 10 times larger than the Thomson pow'er in a 2°C/cm axial 
temperature gradient. 


Table 4.8: Power generated in the melt of GaSb due to 15.7 amp/cm 5 current density. 


Current Density 
amp/cm 2 

Joule 

watt/cm 3 

Thomson 
watt/cm 3 
2°C/cm 20°C/cm 

±15.7 

0.023 

±.00157 ±0.0157 


To measure the thermal perturbations in the vicinity of the liquid-solid inter- 
face. the ampoule was translated into the cooler to promote solidification. The thermocou- 
ple reading was monitored, meanwhile. When the thermocouple reading reached 720 & C the 
translation was terminated. 

This temperature reading was 2-5° above the melting temperature of GaSb. 
i.e. 715°C to 718°C„ With an axial temperature gradient of 20-25°C/cm in the vicinity 
of interface, the thermocouple tip at 720°C temperature reading was expected to be 2 to 3 
mm away from interface into the melt. 

The temperature was measured during application of 10 amp current passed 
for 60 sec from melt to solid and vice versa. Figure 4.65 shows the thermocouple readings 
in the vicinity of the liquid-solid interface, 2 to 3 mm into the melt of GaSb. A 10 amp 
currenl was passed for 60 sec from melt ( + ) to solid (-) for 60 sec. Initially, the thermocouple 
reading decreased from 720°C to 714 f ’C within the first 5-6 sec of pulsation. Afterward, 
the temperature starled increasing gradually. When the pulse was turned off, the cooling 
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effect of the current was terminated and the thermocouple reading started rising to the 
temperature reading before pulsation. A slight overshoot was observed after the current 
was turned off. This could have been due to volumetric Joule heating which had not 
dissipated completely, even after termination of the current. The initial cooling beha\ior 
was due to Peltier and Thomson cooling. Although the cooling effect remained throughout 
the pulse. Joule heating resulted in a gradual rise of temperature. 
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Figure 4.65: Thermocouple readings in a GaSb melt in the vicinity of the liquid-solid inter- 
face during passage of 10 amp current for 60 sec passed from solid, i.e. bottom electrode( + ) 
to melt, i.e. top electrode(-). Before current pulsation, the thermocouple tip was estimated 
to have been 2 to 3 mm away from the interface in the melt. The temperature rise was 
dominated by Peltier heating initially and followed by Joule heating and Thomson heating. 
After termination of current the thermocouple reading returned to the initial steady-state 
value. 

Figure 4.66 shows the thermocouple reading in the melt of GaSb in the vicinity 
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of the liquid-solid interface for a 10 amp current pulse passed for 60 sec from solid(-r) to 
melt (-). These temperature measurements were performed at the same position as the 
profile given in the Figure 4.66, except with different polarity. Initially, the temperature 
started rising rapidly, mostly due to Peltier beating* from 720°C to 726°C within the first 
<5-8 sec of current pulsation and then continued increasing up to 728. 5 D C. After termination 
of current, the temperature decayed to its initial value. 

Figure 4.67 shows temperature measurements similar to those in Figures 4.65 
and 4.66, except the pulse duration was 5 sec. Similar temperature decay and rise was 
observed as in Figures 4.65 and 4.66 during the first 5 sec of the pulsation. 

In these measurements, the transient thermal perturbations in the vicinity of 
the interface was dominated by the Peltier effect and possibly the Thomson effect. Joule 
heating was less pronounced initially. As the current pulsation was continued, the interface 
would have moved to a new position and the Joule heating-induced thermal rise became 
more pronounced. 
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Figu)c 4.66: Thermocouple readings ir. a GaSb melt in the vicinity of the liquid-solid inter- 
face during passage of 10 amp current for 60 sec passed from solid, i.e. bottom electrode(-) 
to melt, i.e. top electrode( + ). Before current pulsation, the thermocouple tip was esti- 
mated to have been 2 to 3 mm away from interface in the melt. The temperature decay 
was dominated by Peltier cooling initially, followed a slight increase due to Joule heating 
and relocation of the interface position. After termination of the current the thermocouple 
reading returned to the initial steady-state value. 
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Figure 4.67: Thermocouple readings in the melt of GaSb charge in the vicinity of liquid- 
solid interface during passage of 10 amp current for 5 sec passed from solid, i.e. bottom 
electrode( + ) to melt, i.e. top electrode(-) and vice versa. Before pulsation, the thermocou- 
ple tip was possibly 2 to 3 mm away from interface into the melt. The temperature rise and 
fall for different polarities were mostly due to Peltier effect. After termination of current 
the thermocouple reading returned to the initial steady-state value. 
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ABSTRACT 

Pb-Sn eutectic alloy was directionally solidified at 1.4 to 3.2 cm/hr with 
forced convection induced by axial vibration of the growth ampoule with a frequency 
of 10 Hz to 40 Hz and an amplitude of 0.5 to 1.0 mm. To determine the exact growth 
rate, an interface demarcation technique was applied. The lamellar spacing was in- 
creased 10% to 40% in ingots solidified with vibration compared to those solidified 
without vibration. The number of grain boundaries was increased by vibration. 

The average intensity of convection in the melt under axial vibration of 
the ampoule was estimated by comparing the experimental results with a theoretical 

model. 


'Current address : State University of Campinas, Campinas, SP, Brazil 

^Current address : MEMC Electronic Materials, 501 Pearl Dr., St. Peters, MO 63376 


PRECEDING -PAGE 


BLANK NOT FILMED 


1. Introduction 

In recent years, many theoretical and experimental investigations have 
been performed on the effect of convection on eutectic microstructure. The main 
motivation for these investigations was that directional solidification of eutectic alloys 
can produce high quality composites with interesting anisotropic properties. 

A theory for diffusion-controlled eutectic growth was presented by Jack- 
son and Hunt [1], Assuming a linear velocity gradient in the flow of melt across the 
interface, Quenisset and Naslain [2], followed later by Baskaran and Wilcox [3] and 
Chandrasekhar et al. [4] calculated numerically the interfacial composition during 
lamellar growth. The variation of composition at the solid/liquid interface was con- 
verted to changes of the undercooling at the interface. Using the extremum criterion, 
a relationship was derived between lamellar spacing and convection. Employing the 
same strategy with two different directions of flow relative to microstructure, Caram 
et al. [5] determined the influence of convection on rod growth. These theoretical 
investigations led to the conclusion that convection increases the lamellar and rod 

spacing during eutectic solidification. 

Eisa and Wilcox [6] and Chandrasekhar [7] experimentally studied so- 
lidification of the MnBi-Bi eutectic with the accelerated crucible rotation technique 
(ACRT). The MnBi rod spacing was increased by using ACRT. 

Popov and Wilcox [8] directionally solidified the Pb-Sn eutectic at 1.0 and 
4.6 cm/hr with ACRT using a vertical Bridgman - Stockbarger technique. Application 
of ACRT did not change the lamellar spacing, although this did influence the rate of 
rotation of the spiral structures observed by Mourer and Verhoeven [9]. 

Another approach to enhance convection in the melt during directional 
solidification is application of vibration. Vibrational mixing has been utilized in melt 
growth by many researchers [10-19], with a wide range of amplitudes and frequencies 
of vibration. For example, ultrasonic vibration at 10 KHz with micron range ampli- 



tudes was used in the Czochralski growth of Te-doped InSb [10-12] and In x Gai_ x Sb 
[13]. The grain size of the Bi 2 Te 3 -Bi 2 Se 3 eutectic system was reduced by ultrasonic 
agitation of the melt during directional solidification [14]. Vibration at frequencies 
of 10 - 100 Hz and an amplitude of 0.05 to 5 cm resulted in better grain selection in 
GaAs [15], CdTe [16-18], and In x Gai_ r Sb [19] crystals directionally solidified by zone 
melting and Bridgman techniques. Even though the vibration enhanced grain selec- 
tion, the number of twin boundaries in CdTe [16-18] and In 0 . 2 Ga 0 .8Sb [19] increased. 
The mechanism for modification of microstructure is not yet known. 

The main objective of this work was to investigate the influence of axial 
vibration of the growth ampoule on the lamellar growth of Pb-Sn eutectic using a 
vertical Bridgman - Stockbarger technique. 

2. Experimental Methods 

The experimental set-up consisted of a vertical Bridgman-Stockbarger 
crystal growth apparatus with provision for axial vibration of the growth ampoule, 
as shown in figure 1. The hot and cold zones were made of Kanthal heating elements 
embedded in Fibrothal insulation, separated by 5 cm of zirconia insulation as an 
adiabatic zone. The ampoule was translated from the hot zone to the cold zone to 
promote solidification. 

The vibration unit consisted of a Bruel Kjaer vibrator connected to a 
power amplifier and HP function generator. With this arrangement the ampoule 
could be oscillated in the axial direction, parallel to gravity, at a frequency of 0.1 
to 100 Hz and an amplitude up to 1.5 mm. The acceleration due to vibration was 
measured using an accelerometer connected to a power amplifier and data acquisition 
system. Figure 2 shows a dynamic acceleration measurement during axial vibration of 
the ampoule at a frequency of 20 Hz and 0.5 mm amplitude. An acceleration of ±0.1 g 
to ±0.15g (g=9.81 m/s 2 ) was measured for the range of frequencies and amplitudes 


utilized for our experiments. 


The growth material was prepared by weighing a proper amount of 99.9999% 
purity Pb and Sn shots, corresponding to 61.1wt% Pb and 38.9wt% Sn. Quartz am- 
poules containing these shots were purged with argon and sealed under a vacuum 
of 10 -6 torr. The materials were alloyed in a rocking furnace at 500°C for 6 hours. 
Then the ampoule was transferred to the vertical Bridgman-Stockbarger furnace for 
directional solidification experiments. Ingots 7 to 12 cm long and 0.6 cm in diameter 
were directionally solidified at an ampoule lowering rate of 1.0 cm/hr rate and an 
axial temperature gradient 40° C/cm (measured using a K-type thermocouple in an 
empty ampoule). The heater and cooler settings were 300°C and 25°C, respectively.- 

Cross sectional samples were taken from several locations along the ingot. 
The samples were cast in a resin mold and mechanically polished. The samples were 
elect rochemically polished in a solution of 800 ml of absolute ethanol, 140 ml of 
distilled water and 60 ml of percloric acid for 60 sec at room temperature. In order 
to reveal the microstructure, the samples were chemically etched in a solution of 1 
part of glycerol, 1 part of acetic acid and 4 parts of nitric acid for 30 sec at room 
temperature. Using such a procedure, the grain boundaries and the lamellar structure 
were revealed. The lamellar structure of cross-sectional and longitudinal samples was 
examined using optical and scanning electron microscopy. 

3. Results and Discussion 

To measure the lamellar spacing, cross-sectional samples were utilized. 
Since the eutectic growth develops with several grains at different orientations, the 
use of longitudinal samples would lead to a mismeasurement of the smaller spacing. 
Figure 3 presents a longitudinal slice. The lamellar spacing appear to be different in 
each grain because of the different orientations relative to the surface. The location 
of a change in the direction of lamella was taken as a grain boundary. 

An intriguing result was obtained from the preliminary experiments per- 
formed without axial vibration. In these experiments, the solidification was carried 


out by moving the ampoule at a constant rate of 1.0 cm/h. The lamellar spacing, 
A, was different from values obtained by others at the same velocity [20-23]. It was 
suspected that the difference between the values of lamellar spacing was due to a 
deviation of the growth rate from the ampoule lowering rate. Sukanek [24] showed 
that the freezing rate may deviate significantly from the lowering rate for directional 
solidification in a Bridgman - Stockbarger furnace, especially near the ends and with 
a large insulation zone (as used here). In our experiments, the sample was not very 
long (ampoule length/ampoule diameter = 12 to 15) and the insulation thickness 
was high (insulation zone/ampoule diameter « 8). Therefore, we concluded that the 
deviation of the freezing rate from the translation rate might have been considerable. 

To assess the deviation of the growth rate from the ampoule lowering 
rate, a technique was used to demarcate the liquid-solid interface periodically during 
the growth period. The growth ampoule was abruptly lowered 0.15 cm every one 
hour. This rapid movement disrupted the structure and enabled the interface shape 
and position to be seen in longitudinal slices. The average macroscopic freezing rate 
was determined by measuring the distance between these interface demarcations and 
knowing the frequency of suddenly lowering the ampoule. Figure 4 shows the growth 
rate versus length fraction solidified in ingots 7.0 to 10 cm long. The freezing rate 
was higher at the first to freeze section of the ingots, decreased halfway through the 
growth, and increased at the end of solidification. These result are in qualitative 

agreement with Sukanek’s prediction [24]. 

Since the growth velocity varied, the lamellar spacing should have also 
varied along the ingots. Figure 5 shows the variation of the lamellar spacing versus 
the fraction solidified. For all vibration conditions, the lamellar spacing was smaller 
near the ends and increased in the center of the ingot, as shown in figure 5. For the 
growth with vibration, the lamellar spacing was larger than without vibration. 

Depending on the intensity of vibration, the use of axial vibration dur- 


